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The Dispersion of Double Refraction in Relation to Crystal 

Structure. 

By T. H. Havelock, M.A., D.Sc., Fellow of St. John's College, Cambridge, 
Lecturer in Applied Mathematics in Armstrong College, Newcastle- 
on-Tyne. 

(Communicated by Professor J. Larmor, Sec. E.S. Received July 5, 1907.) 

In a previous communication,* the writer discussed the double refraction 
of a medium composed of similar particles arranged in rectangular order; the 
results being applied especially to the effects produced in colloidal solutions 
by mechanical stress and by a magnetic field; the present paper is the 
discussion of a similar problem, using a simpler method of procedure 
suggested by Professor Larmor.t The object is to investigate to what extent 
it is possible to consider double refraction, whether produced artificially or 
occurring in natural crystals, as due simply to an aeolotropic distribution of 
similar particles; the questions which arise when one considers whether the 
aaolotropy occurs naturally or is produced ,by the action of mechanical, 
electric, or magnetic force are not specially considered, but The various cases 
are classified in the first section. 

In the next sections we consider a medium composed of a homogeneous 
assemblage of optically isotropic molecules and obtain an expression giving 
the dispersion of the double refraction. This is applied first to artificial 
double refraction produced by mechanical stress and by an electric field; 
then on the same basis the dispersion in quartz is examined. Certain 
divergences in natural crystals and especially such anomalous cases as 
apophylite and vesuvian lead to an extension of the theory. This is, briefly, 
an eeolotropic distribution of optically seolotropic molecules, each of which 
disperses regularly; it is examined first in relation to uniaxal crystals and 
finally in general for biaxal crystals, and is found to contain the possibility 
of all the varieties of dispersion of double refraction which occur in natural 
crystals. 

1. Theories of Double Befr action. 

The ordinary theory of double refraction, in so far as it is molecular, 
amounts to considering the medium as a collection of crystal molecules in 
cubical order, all the varieties of dispersion being postulated of the single 
particle. Thus with the principal axes as co-ordinate axes, if hi) be the 

* ‘Roy. Soc. Proe., 5 A, vol. 77, p. 170, 1906. 

+ Loc. tit., p. 180. 
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polarisation of a typical particle at which the effective electric force is 
(Xi, Y x , Zi), we have equations of the form 

^+7Wi = «Xi, 

-^+h 2 gi = «Yb 

3 J+7'.3%1 = «z x . ( 1 ) 

These equations mean, that the seolotropyis assumed to originate wholly in 
the individual crystal molecule, and that it affects the quasi-elastic force 
under which the polarisation electrons are supposed to vibrate. 

The same method has usually been followed in attempts to apply electron 
theory to the double refraction due to mechanical stress. 

If we consider the equations for an isotropic medium, similar to (1) but 
with a single constant k, then Voigt* assumes that k is altered by amounts pro¬ 
portional to the stresses in the medium; this is a direct effect on the internal 
forces of the molecule, while on the other hand a theory which confines the 
effect to a rearrangement of the molecules in space will express the result 
by a modification of the effective force (Xi, Yi, Z x ) in the previous equations. 

The latter theory appears to comprise the double refraction effects obtained 
by Majorana by placing certain colloidal solutions in a magnetic field. These 
effects appear to be quite distinct from those involved in the Zeemann effect. 
The double refraction obtained in the latter way is explained by a direct 
action of the magnetic field on the vibrating electron, and a result of the 
theory is that the extraordinary and ordinary indices both differ from the 
index of the isotropic medium under no magnetic action; the variation of 
the two indices is in the same direction, one being three times the other, and 
the difference varies inversely as the cube of the wave-length. 

In the Majorana effect, on the other hand, the variations of the two 
indices are in opposite directions, one being approximately twice the other. 

The same distinction occurs in the theory of the Kerr effect, that is, of 
double refraction produced by an electric field. Voigtf has extended his 
theory of the Zeemann effect to cover an analogous effect by an electric field, 
the results being of the same character. Xow Kerr's experiments appeared 
to prove that only the index for vibrations parallel to the field was altered; 
however, in some recent experiments, AeckerleinJ claims to have observed in 

* Yoigt, 4 Annalen der Physik/ vol. 6, p. 461, 1901. 

t Yoigt, ‘Annalen der Pliysik,’ vol. 4, p. 197, 1901. 

j Aebkerlein, 4 Pliysikalische Zeitsclirift, 5 vol. 7, p. 594, 1906 ; also Yoigt, ibid ., p. 811, 
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certain liquids a change in both indices, the variations being in opposite 
senses and one apparently twice the other. It is, of course, possible that 
there may be some direct effect of the electric field upon the periods of the 
electrons, in addition to a merely mechanical effect of rearrangement of 
molecules. 

For the present we are concerned more with working out the possibilities 
of the latter supposition ; that is to say, although Cauchy’s theory of optical 
dispersion, as caused by simple discreteness of the medium, was insufficient 
to account for ordinary dispersion, it is suggested that it may be capable of 
application to the differences of dispersion occurring in different directions in 
a doubly refracting medium. 


2. Optical Properties of a Homogeneous Assemblage of Isotropic Molecules . 

In a dielectric medium, if we write (X, Y, Z) for the electric force, (a, ft, y) 
for the magnetic force, and (/, g , h) for the total electric displacement, we have 
the usual equations 

curl (a, /?, 7 ) = ~ (/, g, h), 

curl(X,Y,Z)=(2) 

The total displacement (f g , h) is given by equations of the form 




/» = Hi/*', (3) 

where ff refers to the polarisation of a typical particle and X* is the number 
of such per unit volume. 

For f h we have an equation 

= “A, (4) 

where X^ is the component of effective electric intensity considered as acting 
on the given particle. 

The determination of the mean value of X^ is the problem in question, and 
in particular, the estimation of the influence of neighbouring particles. 

In the ordinary procedure, we consider round the given particle a sphere 
whose radius B is small compared with the wave-length of the radiation 
crossing the medium, but such that it contains a large number of particles; 
then we evaluate separately the influence of the particles outside and within 
this spherical region. 
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The first part is taken to be the force within a spherical cavity cut in a 
medium uniformly polarised to the value at the point in question of the 
vector 



this, then, contributes a force 


~ {%f h , tg h , th). 


The second part is due to the average effect of the immediately neighbour¬ 
ing particles within the sphere R, and is generally written as 


• 4 tt s 2 tgii, s 3 th), 


where si, s 2 , s 3 are constants depending upon the arrangement of the particles. 
For a simple cubical arrangement, we take the s-constants to be all zero; 
thus, as in the similar argument in the case of liquids, we are able to 
eliminate the consideration of the purely local part of the neighbouring action. 

Now, for any other regular arrangement of particles, the s-constants of this 
method will not all be equal and zero ; but suppose that instead of a sphere R, 
we can choose some other form of surface R', such that we can again ignore 
the local effect of the particles within R', then we are left with the evaluation 
of the force within a cavity R' in a medium uniformly polarised to the value of 
the material polarisation vector at the point in question. In particular, it is 
clear that for a first approximation for a regular arrangement of particles 
differing slightly from a simple cubical system, we may assume R' to be an 
ellipsoid whose axes are nearly equal. 

If we take co-ordinate axes in the directions of the principal axes (a, b, c) 
of the ellipsoid, the force due to the particles outside is now given by 


where 



with similar expressions for B and C. 
Hence, if we put 


(a, b , c) = r (1 + ei, 1 + 6 2 , 1 + e 3 ), 


we have to the first power in 6 X , e 2 , e 3 , 


A = — {1—|6i + f (e 2 + 6 3 )}, 

B = y {1-162 + 1^3 + 6!)}, 

C — ~ { 1 — -f€3+!-(61 + 62)} . 


—7 { 1— l e i + |( e 2 + e 3 )}, 


{I - -t € 3 + |-(61 + 62)} . 


0) 
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Writing these for the moment as 4 tt (± + s) } we have, from (4) and (5), 

^+h% = fx + ^a+s)lf±xjj. ( 8 ) 

If the quantities vary as e ik \ we find from this equation an expression for 
f h ; this can be substituted in equation (3) to give a relation between /and X 
which must necessarily be of the form / = m 2 X/4tt, where n is the refractive 
index. Hence we obtain a value for n which is expressed most conveniently 
in the form 

n 2 — 1 4tt v a oi h n o 2 —1 


— ^ :7r S 

m 2 + 2 + 3s (n 2 —Y) 3 h k h 2 —k 2 


(9) 


^o 2 + 2 ' 

where n 0 is the refractive index of the same substance if isotropic and of the 
same mean density. 

How we substitute from (7) the values for s and expand from (9) to the 
first power in the small quantities ei, e 2 , 63 ; we find the following values for 
the principal refractive indices of the medium :— 

(no 2 -!) 2 


n x — Mo¬ 


ws = no— T V- 


n$ — o%o yj- 


{ 3 ei — (ri-f €2 + 63)}, 
{ 3^2 — (ei + 62 + £3)}* 
{3^3 — (61 + 62+63)}. 


no 

(no 2 -!) 2 
n 0 

(n 0 2 -!) 2 
n 0 

The mean value n 0 is given by 

n 0 = i(n 1 -\-n 2 -!n s ). 

The medium, then, is doubly refracting, and we shall consider first some 
consequences when it behaves like a uniaxal crystal. 

In this case we may take e 2 and e 3 as zero, and we have the extraordinary 
index % and the ordinary index n 2 given by— 

> 0 V- 1 )V 


( 10 ) 


( 11 ) 


ni = w 0 - 


n 0 

i (% 2 —l ) 2 

no 


n 2 = n 0 + e, (12) 

no 

where no = £ (mi 4 - 2n 2 ). 

The mean index, that is the index for the same medium if isotropic and of 
the same density, lies unsymmet rically between the ordinary and extra¬ 
ordinary indices, being twice as far from the latter as from the former. 
Further, if we measure the double refraction S by the difference between the 
two principal indices, we have 

* An 0 2 -!) 2 

o == n 2 —ni = 6 . 

wo 


( 13 ) 

















Downloaded from rspa.royalsocietypublishing.org 


o 
z 
oc 


2<Z}U 


^u 
Sq.o3i/> 



Ss 


w u 

Ko 

H m 



_r a 
< z 

y 5 

b m 
2<zm 

^S=!v5^ 

h^zS 

5 >ui- 

2 i v u 
S0.08 m 



Ss 


w U 

Ko 

H ^ 



1907.] Refraction in Relation to Crystal Structure. 33 

Assuming, then, that the elliptieity of the effective ellipsoid is constant, if 
we know how n 0 varies with the wave-length, this equation gives a law of 
dispersion of the double refraction ; we shall compare this law with some 
experimental data in the following sections. 


3. Double Refraction Due to Mechanical Stress. 


The most suitable test for the above formula would be in the case of 
artificially produced double refraction, as, for instance, when an isotropic 
medium behaves like a uniaxal crystal if subjected to mechanical stress. 
However, in these cases the effect is extremely small, and the principal 
indices ni and n 2 have not usually been measured separately, but simply 
their difference. If the dispersion of the unstrained isotropic medium were 
known either by a curve or a dispersion formula, then it would be possible 
to calculate the dispersion of the difference n\ — n 2 ; for instance, if n 0 , for the 
isotropic medium of the same density, were given by 


then, since 




«+^+ 


c 

\* } 


rs (n 0 2 -- 1) 2 ^ 

6 = U\ — ?i 2 — —-— C, 

no 


we should have S given approximately by a relation of the form 


*-* + S + 5V 


(14) 


Further, the separations of n\ and n 2 from no should be in opposite 
directions, with the former approximately twice the latter. These conse¬ 
quences rest on the initial assumption of the theory, namely, that the 
isotropic medium is, in fact, an assemblage of optically isotropic molecules, 
and not a collection of crystalline molecules, with their axes distributed 
equally in all possible directions. 


4. Electric Double Refraction. 

For the double refraction produced in media by the action of an electric field, 
Kerr observed a dispersive effect, and supposed that the double refraction 
varied inversely as the square root of the wave-length, but this has been 
shown to be erroneous by Blackwell.* His experiments deal with carbon 
disulphide, and curves are given for the variation of the double refraction with 
the wave-length for two different temperatures with a field strength of 
56,000 volts per centimetre. These results give the only data suitable for 
the present purpose. 

* H. L. Blackwell, ‘American Academy Proceedings,’ vol. 41, p. 647, 1906. 
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We have, then, two sets of values for the double refraction 8 of CS 2 at the 
temperatures of 11° C. and 24° C., each set giving values for six different 
wave-lengths with the value of 8 for X = 589 and t = 24° taken as 100 ; 
absolute values can be obtained by means of a determination of 7*04 x 10~ 7 
for 8 at the latter point. 

For the present theory we require the actual values of ni and n 2 in order to 
calculate the mean value n 0 , equal to (ni + 2n 2 )/3 ; however, it is sufficient to 
take for n 0 the refractive index of CS 2 under the same conditions in regard to 
density and temperature. 

To obtain these values for the wave-lengths required, we have a dispersion 
formula for CS 2 at 18° obtained by Martens,* which gives values correct 
to four decimal places at least; the formula is 

o _ , Tfb K 2 7 a 2 

^ = m+ xCT~ /cX ’ 

where m = 1*63403; h = 0*000300; 

m'= 0*88113; X' = 217*53. 

The correction for temperature can be made to sufficient accuracy by 
allowing a decrease of index of 9 x 10~ 4 for each degree rise of temperature.! 
Making these calculations, we have the values of n 0 used in the following 
Tables I and II. The first column gives the wave-length X, the second n 0 , 
the third the observed values of the double refraction 8; the fourth column 
gives the calculated values of S?i 0 /(?i 0 3 — l) 2 , which, according to the theory, 
ought to be constant. 


Table I.—Carbon Disulphide at 24°. 


A. 

*0- 

<^obs. = n i —w 2 . 

C = 5» 0 /( w 0 a “ l) 2 - 

^calc. 

421 *6 

1*678 

119 *4 

60 *78 

118 *8 

441 *5 

1 *668 

115 *0 

| 60 *39 

115 *2 

467 *7 

1*656 

110 *5 

60 *27 

110 *9 

508 *5 

1 *642 

106 *0 

60 *50 

106 *0 

538 *7 

1 *634 

103 *0 

60 *35 

103 *2 

589 *0 

1*624 

100 *0 

60 *58 

99 *8 


Mean value of 0 = 60*48 ; greatest divergence = 0*5 per cent, approx. 


The fifth column in each table was obtained by using the mean value of 
0, together with the values of i n order to recalculate the values of 8. A 
comparison of the third and fifth columns shows the amount of agreement 
between the observed and calculated values of 8; it may be remarked that 

* Martens, ( Annalen der Physik,’ vol. 6, p. 632, 1901. 
t Landolt, ‘Tabellen, 5 3rd edition, p. 671. 
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1907.] Refraction in Relation to Crystal Structure . 35 


the observed values are the mean values of results varjung by as much as 
2 per cent. 

Table II.-—'Carbon Disulphide at 11°. 


A. 

n 0 . 

<5obs. = n l — n 2 . 

C = Sn 0 j(n 0 2 —1) 2 . 

$calc. 

421 '6 

1 -689 

125 '6 

61 '80 

125 '7 

441 '5 

1 -679 

121 -0 

61 *40 

121 *8 

467*7 

1-667 

116 *8 

61 *62 

117 *2 

508 *5 

1 -653 

112 *4 

62 '01 

112 '1 

538 '7 

1 -645 

109 *4 

61 *93 

109 '2 

589 '0 

1 -635 

106 *6 

62 *26 

105 *9 


Mean value of C = 61'84 ; greatest divergence = 0*7 per cent, approx. 


We can calculate the absolute value of C from the value for 8 at t = 24° 
and A = 589, viz., 7*04 x 10“" 7 . If, then, e is the ellipticity of the effective 
ellipsoidal cavity of the theory, we have e equal to 50, and we find the 
following values:— 

e = 0*000,002,18 at 11°, 

e = 0*000,002,13 at 24°. (15) 


5. Dispersion of Double Refraction in Quartz. 


In order to determine to what extent the double refraction of natural 
crystals can be accounted for by an aeolotropic distribution of isotropic 
molecules, we require data for some simple crystal for which the difference of 
the indices is sufficiently small; for this purpose we have accurate determina¬ 
tions of and n 2 for quartz extending over a very wide range of wave-lengths. 
The values for these given in Table III are by Martens* for wave-lengths 
ranging from 185 to 768, and by Carvallof from 801 to 1849. 

The numbers in the fifth column C are calculated in each case from 
the relations 


8n 0 

(V-l) 2 ’ 


n 0 = \(ni + 2n 2 ). 


Then the quantity C ought to be a constant. Kemembering that we have 
neglected terms of the second order, and omitting for the present the extreme 
limits in the ultra-violet and infra-red, we see that the figures give 
considerable support to the theory. For throughout the large range of X, 
from 358 to 1159, C does not differ from a mean value of 0*007245 by as 
much as 0*2 per cent, in any case. For the sake of comparison one might 

* Martens, ‘ Annalen der Physik,’ vol. 6, p. 628, 1901. 
t Carvallo, ‘ Comptes Rendus,’ vol. 126, p. 730, 1898. 

VOL. LXXX.—A. 
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calculate the values of some other relation which might suggest itself; for 
instance, the quantity Sno/(n 0 2 — 1) decreases continuously throughout the 
whole range, its extreme values differing from a mean value by as much as 

16 per cent. 

Table III.—Quartz. 


X. 

%• 

» 2 . 

1 

0 X 10°. 

185 

1 *68988 

1 -67571 

0 *01417 

7151 

198 

6394 

5087 

1307 

7150 

224 

3231 

2045 

1186 

7175 

242 

1647 

0522 

1125 

7170 

263 

0386 

1*59306 

1080 

7174 

281 

1 *59516 

8467 

1049 

7188 

308 

8564 

7547 

1017 

7204 

340 

7737 

6747 

0 -00990 

7223 

358 

7369 

6390 

979 

7236 

394 

6805 

1 -55846 

959 

7235 

434 

6339 

396 

943 ■ 

7235 

467 

6034 

102 

932 

7231 

486 

1 -55897 

1 *54967 

930 

7253 

508 

746 

822 

924 

7246 

533 

599 

680 

919 

7246 

546 

534 

617 

917 

7245 

560 

461 

546 

915 

7251 

589 

335 

424 

911 

7252 

627 

186 

280 

906 

7253 

643 

131 

227 

904 

7253 

656 

091 

189 

902 

7247 

670 

047 

145 

902 

7258 

768 

1 -54794 

1*53903 

891 

7238 

801 

725 

834 

891 

7257 

832 

661 

773 

888 

7250 

867 

598 

712 

886 

7250 

905 

532 

649 

883 

7245 

946 

464 

583 

881 

7245 

991 

392 

514 

878 

7240 

1042 

317 

442 

875 

7246 

1097 

238 

366 

872 

7232 

1159 

152 

283 

869 

7232 

1229 

057 

192 

865 

7224 

| 1307 

1 -53951 

090 

861 

7218 

1396 

832 

1*52977 

855 

7202 

1 1497 

692 

842 

850 

7198 

1 1615 

524 

679 

845 

7201 

| 1749 

319 

485 

834 

7164 

i 1849 

163 

335 

828 

7154 


In consequence of the use of Babineffs compensator in determining the 
double refraction of minerals, the dispersion of double refraction in quartz has 
been studied extensively. From experimental observations in the visible and 
ultra-violet spectrum, Mace de Lepinay* found that S could be represented very 
accurately by a formula— 


^ . b , c 


* Mace de L6pinay, ‘ Annales de la Faculte des Sciences de Marseille,’ vol. 50, p. 1, 
1891. 
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The numerical values of the constants were found by Carvallo to fit also 
his observations in the infra-red. This relation is of the form to be expected 
from our formula connecting S and n 0 , assuming no to be expressible by a 
dispersion formula of the usual type. 

We have assumed in the calculations that the effective cavity was an 
•ellipsoid of small ellipticity. We might then examine the order of magnitude 
of the second order terms which have been neglected, but it is simpler to 
obtain now a more general relation free from this approximation. 

Using the notation of equations (5) to (9), the only assumption we make 
about the effective cavity is that the force within it, due to the polarised 
medium outside, has components 

^K 1 )' 

z), (16) 

where Si, s 2) s 3 are constants. 

Then assuming as before that the molecules are isotropic, we have a 
principal index given by 

_ n 2 — 1 ___ n 2 — 1 

^ + 2 + 3^ 2 -l) “ no 2 + 2 * 

The quantity (n 2 -f- 2)/(n 2 —1 ) differs by a constant from its value in any 
other state; or, more simply, we have 


1 _ 1 

n 2 — 1 no 2 — 1 

Consequently, if the medium is uniaxial, we have the relation 

——— ——= constant, 
n 2 2 — 1 ni 2 — 1 


(17) 

(18) 


where ni and n 2 are the extraordinary and ordinary indices. 

Using the values given in Table III, we can examine this relation for 
quartz ; the result is shown in Table IY. 


Table IY.—Quartz. 


A. 

C 

O'. 

A. 

358 

0 *01439 

0 *01443 

656 

434 

0 *01439 

0 *01439 

768 

486 

0 *01442 

0 *01442 

832 

533 

0 *01441 

0 *01442 

905 

560 

0 *01442 

0 *01440 

991 

627 

0 *01442 

0 *01439 

1097 


O' = (n 2 2 — l) -1 — (n^ — 1) 1 . Mean value = 0 *01441. 
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6 . Some Exceptional Crystals . 

The relation obtained in (18) does not cover all the cases of regular dis¬ 
persion, as we see from the figures for Iceland Spar given in Table V. But 
it may be noticed that in this case the double refraction 8 is 20 times 
greater than in quartz ; thus simple distribution of isotropic molecules is 
apparently not quite sufficient to account for so large a difference between 
the two indices. 


Table Y.—Iceland Spar. 


A. 

II 

c'^cv-i)" 1 -^ 2 -!)- 1 . 

410 

0 *1837 

0 *2581 

434 

0 -1812 

0 -2578 

508 

0 -1757 

0 -2565 

589 

0 -1719 

0-2555 

656 

0-1698 

0 ‘2548 

768 

0 *1671 

0 -2538 


Further, we have exceptional cases in which 8 decreases with the wave¬ 
length, but much more rapidly than in other crystals. For instance, in 
strontium hyposulphate, the rate of decrease is about 10 times that of 
quartz, and we find that the quantity C of the previous calculations is by no- 
means constant, but decreases with the wave-length. 

If we have 8, given by the relation 


we see that the substance may be either a positive or negative uniaxal crystal,, 
according as e is negative or positive ; but in either case the absolute value 
of 8 decreases with n 0 , that is, in general, 8 will decrease as X increases. How 
this is the rule in most actual crystals, but there are a few for which 8 
increases with X, although no diminishes; the figures for apophylite 
are given in Table VI for the purpose of comparison with a later set of. 
figures. 

Table VI.—Apophylite.* 


A. 

8. 

486 -1 

0 -00218 

501 *6 

0 -00232 

546 -1 

0 *00241 

587 -6 

0 -00265 

589 -6 

0 -00265 

656 -3 

0 -00284 

667 -8 

0 -0C)289 


* B. Trolle, ‘ Physikalische Zeitschrift,’ vol. 7, p. 700, 1906. 
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1907,] Refraction in Relation to Crystal Structure . 39 

Also, there are crystals for which h passes through a minimum, and some 
for which the double refraction changes from positive to negative—the sub¬ 
stance being isotropic for some definite wave-length. 

It appears, then, that a theory of seolotropic distribution of isotropic mole¬ 
cules which may serve for the artificial double refraction produced in isotropic 
bodies and for simple crystals is not sufficient to cover the varieties of dis¬ 
persion found in natural crystals. This is naturally to be expected; for, in 
addition to the varieties of dispersion in biaxal crystals, to be considered 
later, we see that if such a simple theory were sufficient, then the symmetry 
of all the physical properties of a crystal would be the same as its optical 
symmetry, and this is not the case. The direction in which to modify the 
theory is clear; we shall consider a crystalline medium as a homogeneous 
assemblage of crystal molecules. Whether or no the crystal molecule can 
be identified with the chemical molecule need not be considered here. 
Leaving the general theory till later, we shall consider first the case of a 
uniaxal crystalline medium. 


7. Modified Theory of Uniaxal Crystals . 

The previous theory can be extended in the following manner:—The 
molecules (or crystal units) of the medium are not necessarily to be supposed 
ellipsoidal in shape, but are optically aeolotropic, so that the subsidiary 
equations connecting the polarisation of a particle with the effective electric 
force are seolotropic, with an axis of symmetry; the particles are supposed to 
be arranged in a homogeneous assemblage such that the effective cavity may 
be taken as an ellipsoid of revolution of ellipticity e, and having its axis 
of symmetry coincident in direction with that of the crystal unit. Then, 
instead of equation (9), we have now 

ni 2 —1 _ 47r y g h __ ??, 01 2 -— 1 

% 2 + 2 + 3s 1 (w 1 2 -l) 3 /l h' 2 -k 2 W 2 + 2’ V ' 


with a similar equation in and n 02 . 

In these, n 0 1 and n 0 2 are the extraordinary and ordinary indices for the 
substance, with a regular cubical arrangement of the same molecules in the 
same density. 

Hence, substituting for Si and their values in terms of e, we have the 
two principal indices given by 


ni = ?ioi 


2 

15 


OV-I ) 2 

•^01 


e, 


n 2 = %02 + tV 


(no2 2 


-i ) 2 


^02 


( 20 ) 
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And the double refraction 8 is given by 

* ' i f2(noi 2 -l) 2 fn 02 2 ~-l) 2 \ 

o == Woi — ^02 “tV 1 ——--+ --- [ 6. 

L ^01 ^02 J 

These equations (20) and (21) would represent the effect of deformation in 
a natural uniaxal crystal. For the present purpose we apply them to a 
natural unstrained crystal. They cannot be verified by actual calculation in 
this case, but we shall see that they cover all the anomalous cases mentioned 
in the previous section. 

To find how 8 varies with the wave-length, we have, from (21), 


dS_ 

f I 2 ( 3m 01 2 + 1) Ool 2 — 1) / 

o 

>-• 

1 

, , (3wo2 2 +1)Oo2 2 — 1)/ 

[ 3«02 

dX 

1 ^ 

f 3x L 

15 U02 2 

idx' 


( 22 ) 


1STow the object is to account for anomalous dispersion of 8 in regions free 
from absorption. Then if n 0 1 and n 02 follow the regular law of diminishing 
as X increases, we have both dn 01 /d\ and dn^/dx negative; and as we can 
have n 0 i greater or less than n 02 and e positive or negative, we see that it is 
possible to have 3S/0X positive. In fact, it is possible to have the following 
varieties of change of 8 as X increases : (i) slow decrease, (ii) rapid decrease, 
(iii) slow increase, (iv) S passing through a minimum, and (v) 8 changing in 
sign. 

These cases all occur, being typified by quartz, strontium hyposulphate, 
apophylite, and for the latter cases by varieties of vesuvian and melilite.* 

The theory can be illustrated best by building up from a known crystal of 
regular dispersion of double refraction one which behaves like apophylite. 
Suppose we consider the uniaxal crystal phenakite, for which there are 
available determinations of the principal indices for five different wave¬ 
lengths; we denote these indices by n 0 i and n 0 2 , and we find that their 
difference diminishes as the wave-length increases. Suppose, now, that the 
crystal molecules of phenakite are rearranged so that the effective cavity has 
an ellipticity of 0*035, and has its axis along the direction of the axis of 
symmetry of the crystal units; then by equation ( 21 ) we may calculate the 
double refraction 8 of the medium thus built up. The results are shown in 
Table VII. We see that although the double refraction of the phenakite 
crystals follows the regular law of decreasing as X increases, yet the resulting 
double refraction 3 of the distorted medium increases as X increases, and is of 
about the same order as that of apophylite (Table VI). This is, of course, 
only intended as an illustration of the possibility of explaining the irregular 
dispersion in regions which are not in the vicinity of absorption bands. 

* Hlawatseh, ‘Tschermak’s Mitteilungen,’ vol. 23, p. 4X5, 1904. 


[July 5, 

( 21 ) 
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Table VII. —Phenakite with e = 0*035. 
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508 

537 

589 

643 

670 


1 *67451 
1*67254 
1 *66977 
1 *66735 
1*66639 


1 *65858 
1 *65664 
1*65394 
1*65154 
1*65060 


0 *01593 
0 *01590 
0 *01583 
0 *01581 
0 *01579 


0 *00255 
0 *00260 
0 *00265 
0 *00273 
0*00275 


8 . Dispersion in Bicixal Crystals . 


We consider now the previous theory in the most general form, namely* 


the double refraction of a medium composed of optically seolotropic particles 


with the same orientation and arranged so as to give an ellipsoidal effective 
cavity with its principal axes in any three fixed mutually perpendicular 
directions; we shall indicate briefly the general character of the results. 

We take the principal axes of a typical particle as co-ordinate axes, and 
have (\i, n i, vi) y (\ 2 , /x 2 , v 2 ), and (Vs, ^ 3 , z> 3 ), as the direction cosines of the 
principal axes of the effective cavity. Then if (f, g\ hf) is the polarisation 
of the medium, the ^-component of the force within the cavity due to the 
polarisation is given by 


Xi = (\i 2 A 4 X 2 2 B 4* X 3 2 C)/ / 4 (X 1 M 1 A 4 X 2 / 4 2 B 4 Xg/^C) g f 


4 (\ 1 V 1 A 4 X 2 z/ 2 B 4 X 3 v 3 0) h\ (23} 


where A, B, C, are the usual constants for the cavity and can be written as— 


(A, B, C) = — 4 4?r ($ 1 , s 2 , S 3 ). 


Further, if we write n 0 1 , n 02 , n 0 3 for the principal indices of a medium 
composed of the same units, but arranged in regular cubical order, we have 
three equations of the form 


^ 2 1 a _ 



(24) 


We can obtain now three equations connecting the components of electric 
force (X, Y, Z) with the polarisation (/', g* h '); but if (/', m f , n') are the 
direction cosines of a principal axis of the resulting ellipsoid of polarisation 
and if n is the corresponding principal index of refraction for this direction, 
we must have 


Vf + m’c/ + n’h' = (I'X + m'Y + n'Z). 


(25) 


Hence, for the determination of the principal axes and indices of the 
medium we obtain three equations of the type 



l t j r cii 2 m f 4 ctitfi = 0, 


(26) 
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where the a-constants are given by relations of the form 

ail = \l 2 Si + ^2^2 + ^3 2 ^3 j ^23 = /^1 Vl^I +/^2^2S2 +/^3^3 S 3- (27) 


We shall examine now the various types of crystals covered by these 
equations. 

(a) Prismatic system .—If we take the simplest case, we have the principal 
axes of the crystal unit coinciding with those of the effective cavity; then 
-evidently the three principal polarisation axes are in the same directions, and 
maintain these directions for all wave-lengths, as in crystals of the prismatic 
type. 

From equations (26) and (27) we see that the principal refractive indices 
are given by 


Thus we have 



(29) 


Hence, as in the theory of § 7 for uniaxal crystals, we may have the 
indices varying differently with the wave-length and consequently the 
possibility of such crystals as brookite, in which the order of magnitude of the 
indices changes and the plane of the optic axes is turned through a right angle. 

(b) Monoclinic System. —Suppose, now, that one principal axis of a 
particle coincides with a principal axis of the cavity, then this will be an 
axis of the polarisation ellipsoid which is fixed in direction; thus we obtain 
the monoclinic system of crystals with a single plane of symmetry. 
Without discussing the equations in detail, we consider, as an illustration, the 
dispersion of the other two polarisation axes. These are at right angles to 
each other and lie in the plane of symmetry of the crystal. Let <£ be the 
angle giving their position with reference to the fixed co-ordinate axes, and 
a the corresponding angle for the position of the fixed axes of the effective 
cavity in the same plane. 

Then, supposing the plane of symmetry to be the (2, 3) plane, we find, 
after substitution in (26) and (27) and reduction, that 


(^ 2 —S 3 ) sin 2ol 


tan 2 (f) 


(30) 


(s 2 — s 3 ) cos 2a— (n 02 2 — 1) 1 + (% 3 2 — 1) 


Now s 2 , S 3 , and are assumed to be independent of the wave-length; suppose, 
then, that the indices n 02 and %>3 for the crystal molecules in cubical order 
decrease regularly with the wave-length, so that on the simplest supposition 
we have dispersion formulae of the type 


2 _ 


b 
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1907.] Refraction in Relation to Crystal Structure. 43 


Then from (30) we see that the dispersion of the angle 0 is given by 
a relation of the form 


tan 20 = 


p + gX 2 + r \ 4 
V + ^X 2 + ? y X 4 ' 


(31) 


Sufficient experimental data on the dispersion of the principal axes of 
monoclinic crystals are not available to test this formula; but a similar 
relation without terms in X 4 has been suggested on other grounds by 
Nakamura,* and found to agree fairly well with some cases. 

(c) Anorthic System .—Finally we have the general case in which none of 
the axes of the crystal unit coincides with a principal axis of the effective 
•cavity, and consequently the three principal axes of polarisation all change, 
both in magnitude and direction, with the wave-length; this includes, then, all 
.anorthic crystals. 


9. Summary. 

Beginning with a theory intended to apply to artificial double refraction 
produced in an isotropic medium, we found that the consequences could be 
expressed simply in the following manner:—If n is a refractive index of the 
substance, the quantity is altered by an amount independent of the 

wave-length, and varying only with the direction; this implies a small 
change in n which is proportional to (;x 0 2 —l)/^o, where n 0 refers to a standard 
state of arrangement in regular cubical order. 

Passing to natural crystals, we find anomalous dispersion of the double 
refraction occurring in regions away from absorption bands, and a theory 
modified to account for these cases leads to a consideration of the optical 
relations of crystal structure. In this connection it is of interest to compare 
some recent work in two directions. On the one hand we have attempts to 
explain crystal structure as a function of chemical constitution; the crystal 
is described as a homogeneous assemblage of spheres which are combined into 
sub-groups defining the crystal unit or molecule. On the other hand, we 
have work connecting the dispersion of double refraction in certain groups of 
minerals with their chemical constitution, some members being found to 
possess a minimum of double refraction, f 

Consequently, it should be possible to explain the varieties of dispersion 
of double refraction by means of the structure of the crystal, and we have 
attempted this in the previous sections, with the following assumptions:— 
The crystal unit contains vibrating electrons, so that their combined effect is 
expressed by three principal equations connecting the polarisation of the 


* Nakamura, ‘ Physikalische Zeitschrift,’ vol. 6, p. 172, 1905. 
t Hlawatsch, loc. cit. 
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4 4 Dispersion of Double Refraction in Relation to Crystal Structure. 

unit with the effective electric field; if, then, these units are arranged in 
regular cubical order, we have a medium with principal refractive indices 
along three fixed directions in space, and we assume that in this case 
we have regular decrease of the double refraction with increasing wave¬ 
length in regions away from absorption bands. But considering in general 
any other homogeneous assemblage, we express the effect by a change in the 
effective electric field acting on the crystal unit; this effect we have 
estimated by supposing, as a sufficient approximation, that the effective 
cavity is slightly ellipsoidal instead of being spherical. Thus differences of 
packing of the crystal molecules are represented optically by variations in 
the ratios of the axes of the effective cavity and in their directions in space- 
compared with the polarisation axes of the individual unit. Combining 
these assumptions, we have shown that they are sufficient for a descriptive- 
theory covering the varieties of dispersion of double refraction found in 
natural crystals. 
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